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Abstract 



The Yukawa couplings of the tau lepton and the bottom quark become 
comparable to, or even exceed, electroweak gauge couplings for large values 
of the SUSY parameter tan /3. As a result, the lightest tau slepton fi and 
bottom squark bi can be significantly lighter than corresponding sleptons and 
squarks of the first two generations. Gluino, chargino and neutralino decays to 
third generation particles are significantly enhanced when tan /3 is large. This 
affects projections for collider experiment reach for supersymmetric particles. 
In this paper, we evaluate the reach of the Fermilab Tevatron pp collider 
for supersymmetric signals in the framework of the mSUGRA model. We 
find that the reach via signatures with multiple isolated leptons (e and fj.) is 
considerably reduced. For very large tan /3, the greatest reach is attained in 
the multi-jet+-9y signature. Some significant extra regions may be probed by 
requiring the presence of an identified 6-jet in jets+^^n events, or by requiring 
one of the identified leptons in clean trilepton events to actually be a hadronic 
1 or 3 charged prong tau. In an appendix, we present formulae for chargino, 
neutralino and gluino three body decays which are valid at large tan /?. 
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Typeset using REVT^ 



I. INTRODUCTION AND MOTIVATION 



The minimal supergravity model (mSUGRA) is commonly regarded as the paradigm 
framework for phenomenological analyses of weak scale supersymmetry. The visible sector 
is taken to consist of the particles of the Minimal Supersymmetric Standard Model 
(MSSM). One posits, in addition, the existence of "hidden sector" field(s), which couple to 
ordinary matter fields and their superpartners only via gravity. The conservation of i?-parity 
is assumed. Supersymmetry is broken in a hidden sector of the theory; supersymmetry 
breaking is then communicated to the visible sector via gravitational interactions. The 
technical assumption of minimality implies that kinetic terms for matter fields take the 
canonical form; this assumption, which is equivalent to assuming an approximate global 
U{n) symmetry between n chiral multiplets, leads to a common mass squared ttiq for all 
scalar fields, and a common trilinear term Aq for all A parameters. These parameters, 
which determine the sparticle-particle mass splitting in the observable sector are taken to 
be comparable to the weak scale, M^^ak ■ In addition, motivated by the apparently successful 
gauge coupling unification in the MSSM, one usually adopts a common value mi/2 for all 
gaugino masses at the scale Mqut — 2 x 10^^ GeV. For simplicity, it is commonly assumed 
that in fact the scalar masses and trilinear terms unify at Mqut as well. The resulting 
effective theory, vahd at energy scales E < Mqut, is then just the MSSM with the usual 
soft SUSY breaking terms, which in this case are unified at Mqut- The soft SUSY breaking 
scalar and gaugino mases, the trilinear A terms and in addition a bilinear soft term B, 
the gauge and Yukawa couplings and the supersymmetric yU term are all then evolved from 
Mqut to some scale M ~ M^eak using renormalization group equations (RGB's). The large 
top quark Yukawa coupling causes the squared mass of one of the Higgs fields to be driven 
negative, resulting in the breakdown of electroweak symmetry; this determines the value 
of /i^. Finally, it is customary to trade the parameter B for tan/5, the ratio of Higgs field 
vacuum expectation values. The resulting weak scale spectrum of superpartners and their 
couplings can thus be derived in terms of four continuous plus one discrete parameters 

'^o, '"^1/2, ^0, tan/3 andsgn(/i), (1.1) 

in addition to the usual parameters of the standard model. 

The consequences of the mSUGRA model have been investigated for collider experiments 
at the CERN LEP2 e+e" collider g, the Fermilab Tevatron pp collider [|,|, the CERN 
LHC pp collider and a possible next hnear e+e~ collider (NLC) operating at -^i ~ 500 
GeV . In all but the last of these studies (where the effect of the tau Yukawa coupling on 
aspects of the phenomenology of the stau sector is carefully examined), small to moderate 
values of the parameter tan/? ~ 2 — 10 have been adopted. This was due in part to the 
fact that event generators such as ISAJET had not been constructed to provide reliable 
calculations for large tan (3. In particular, effects of tau and bottom Yukawa couplings, 

V2Mwcosl3' V2Mwcos(3 ^ ' ' 

which become comparable to the electroweak gauge couplings and even to the top Yukawa 
coupling ft = grrit/ {V2M w sin P) if tan/3 is large, had not been completely included. The 
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correct inclusion of these couplings has a significant impact [1T0| , [TT| on the search for super- 
symmetry at colliders. 

In the mSUGRA model, the parameter tan/3 can be as large as tan/? ~ mt/mh, where 
the quark masses are evaluated at a scale ~ M^eak'-, since the running m;, is considerably 
smaller than 5 GeV, tan (3 values up to 45-50 are possible. Such large tan (3 values are indeed 
preferred in some 5*0(10) GUT models with Yukawa coupling unification. In practice, one 
finds that if tan/3 is chosen to be too large, fb diverges before Mqut- A slightly stronger 
upper limit on tan (3 is obtained from the requirement that the mass of the pseudo-scalar 
Higgs boson, should be positive. The precise value of the upper bound on tan /? depends 
somewhat on the other mSUGRA parameters. 

In a recent Letter fill], we reported on an upgrade of the event generator ISA JET that 
correctly incorporated the effects of r and h Yukawa interactions so that it would provide 
reliable predictions for supersymmetry with large tan/3. Novel phenomenological implica- 
tions special to large values of tan/? were pointed out: in particular, it was noted that while 
Tevatron signals in multilepton (e and /i) channels were greatly reduced, there could be new 
signals involving 6-jets and r-leptons via which to search for SUSY. In this paper, we focus 
our attention on the search for supersymmetry at the Main Injector (MI) upgrade of the 
Fermilab Tevatron pp collider, (a/s = 2 TeV, integrated luminosity / Cdt = 2 fb~^) and the 
proposed TeV33 upgrade (a/s = 2 TeV, integrated luminosity / Cdt = 25 fb~^) for the case 
where tan /3 is large. 



A. Sparticles masses at large tan/3 

Large b and r Yukawa couplings significantly alter the mass spectra of the sparticles and 
Higgs bosons as shown in Fig. 1. Here we plot various sparticle and Higgs boson masses 
versus tan/3 for mSUGRA parameters mi/2 = 150 GeV, Aq = and a) itiq = 150 GeV and 
b) rriQ = 500 GeV, for both signs of /i. We fix the pole mass = 170 GeV. 

The b and r Yukawa couplings contribute negatively to the renormalization group running 
of the sbottom and stau soft masses, driving them to lower values than soft masses for 
the corresponding first and second generation squarks and sleptons. In addition, the off- 
diagonal terms in the sbottom and stau mass-squared matrices mb{—Ab + /xtan/3) and 
mr{—Ar + /itan/3) can result in significant mixing between left and right sbottom and stau 
gauge eigenstates, and a possible further decrease in the physical masses for the lighter of 
the two sbottom (and stau) mass eigenstates m^^ and rrif^. If tan/3 is small, fi ~ f^, while 
(because of top quark Yukawa interactions) bi ^ bi- The impact of bottom and tau Yukawa 
interactions can be seen in Fig. 1: ruf-^ ~ rrien ^"^^ tan/3, and as tan/3 increases, m^^ 
decreases, while m^^ remains constant. Likewise, decreases with increasing tan/3, while 
mj^ remains constant. In the case of frame a), ultimately m^^ drops below Tny^^^ and m^^ 

so that the two body decays Wi fiUr and Z2 tit become allowed, and dominate the 
branching fractions. 

It is well known that at low to moderate values of tan /3, the large top Yukawa coupling 
drives the Higgs mass m^^^ negative values, resulting in a breakdown of electroweak 
symmetry. At large tan/3, the large b and r Yukawa couplings drive the other soft Higgs 
mass-squared m|^^ to small or negative values as well. This results overall in a decrease in 
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mass for the pseudo-scalar Higgs tua relative to its value at small tan/3. Since the values 
of the heavy scalar and charged Higgs boson masses are related to tua, they decrease as 
well. This effect is also illustrated in Fig. 1, where the mass decreases dramatically with 
increasing tan /3. The curves are terminated at the value of tan f3 beyond which m\ < 0, and 
the correct pattern of electroweak symmetry breaking is not obtained as already mentioned. 
We found that the pseudoscalar mass niA, obtained using the 1-loop effective potential, is 
unstable by up to factors of two against scale variations for relatively low values of scale 
choice Q ~ Mz- This instability would be presumably corrected by inclusion of 2-loop 
corrections. We find the choice of scale Q ~ ^JfrTiJrTi^ to empirically yield stable predictions 
of Higgs boson masses in the RG improved 1-loop effective potential (where we include 
contributions from all third generation particles and sparticles). This scale choice effectively 
includes some important two loop effects, and yields predictions for light scalar Higgs boson 
masses in close accord with the results of Ref . |]l2l . 



B. Sparticle decays at large tan (3 

For large values of tan/5, h and r Yukawa couplings become comparable in strength to 
the usual gauge interactions, so that Yukawa interaction contributions to sparticle decay 
rates are non-negligible and can even dominate. This could manifest itself as lepton non- 
universality in SUSY events. Also, because of the reduction of masses referred to above, 
chargino and neutralino decays to stau, sbottom and various Higgs bosons may be allowed, 
even if the corresponding decays would be kinematically forbidden for small tan/? values. 
The reduced stau, sbottom, and Higgs masses can also increase sparticle branching ratios 
to third generation particles via virtual effects. These enhanced decays to third generation 
particles can radically alter the expected SUSY signatures at colliders. 

We have re-calculated the branching fractions for the g, hi, ii, fj, u^-, Wi, Zi, h, H, A 
and if^ particles and sparticles including sbottom and stau mixing as well as effects of 



b and r Yukawa interactions. For Higgs boson decays, we use the formulae in Ref. |]13 
We have recalculated the decay widths for g tbWi and g hhZi. These have been 
calculated previously by Bartl et al. ; our results agree with theirs if we use pole fermion 
masses to calculate the Yukawa couplings. In ISAJET, we use the running Yukawa couplings 
evaluated at the scale Q = nig (nit) to compute decay rates for the gluino {Wi,Zi). This 
seems a more appropriate choice, and it significantly alters the decay widths when effects of 
fb are important. The Zi —>■ rfZj and Zi —>■ bbZj decays take place via eight diagrams (/i,2, 
2, Z, h, H and A exchanges). In our calculation of g and Z^ decays, we have neglected b 
and T masses except in the Yukawa couplings and in the phase space integration. We have 
also computed the widths for decays Wi ZjTu which are mediated by W, fi^2, t^t and 
exchanges; in these cases, we retain effects only in the Yukawa couplings. Formulae for 
these three-body decays are presented in the Appendix. 

To illustrate the importance of the Yukawa coupling effects, we show selected branching 
ratios of Wi and Z2 in Fig. 2. In all frames we take /i > 0. Frames a) and h) are for 
the mSUGRA case (mo, mi/2, Aq) = (150,150,0) GeV; frames c) and d) show the same 
branching fractions, but take mo = 500 GeV instead. In frame a), for low tan/3 we see that 
the Wi euZi and Wi — > TuZi branching ratios are very close in magnitude, reflecting the 
smallness of For tan/5 ^ 10, these branchings begin to diverge, with the branching to 
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r's becoming increasingly dominant. For tan (3 > 40, the two body mode Wi — ^ fiu opens 
up and quickly dominates. Since this decay is followed by fi — > tZi, the end products of 
chargino decays here are almost exclusively tau leptons plus missing energy. 

In frame b), we see at low tan/3 the Z2 eeZi and Z2 — > rfZi branchings are large 
(~ 10%) and equal, again because of the smallness of the Yukawa coupling. Except for 
parameter regions where the leptonic decays of Z2 are strongly suppressed, W1Z2 production 
leads to the well known 3i (= e, fi) signature for the Tevatron collider |T5| . As tan /? increases 
beyond about 5, these branchings again diverge, and increasingly Z2 TfZ\ dominates. 
Results of phenomenological analyses of trilepton signals for tan /5 ~ 8 — 10 obtained using 
older versions of ISAJET should, therefore, be interpreted with caution. For tan /3 > 40, 
Z2 — > rfi opens up, and becomes quickly close to 100%. Near the edge of parameter space 
(tan/5 ~ 45), the Z2 — > Zih decay opens up, resulting in a reduction of the Z2 — > rfi 
branching fraction. 

In frame c), the large value of mg = 500 GeV yields a large value of nif^ (and other slepton 
masses) even if tan /3 is large. In this case, the W\ branching fractions are dominated by the 
virtual W boson, so that B{yV\ Zxtv) and B(yV\ ZiTu) are nearly equal over almost 
the entire range of tan/5. The branching fractions of Z2 for mo = 500 GeV are shown in 
frame d). As in frame c), the branching fraction of Z2 to r's and e's is nearly the same 
except when tan/? > 35 — 40. In this case, there is a steadily increasing branching fraction 
of Z2 — > Zibb (and to some extent, also of Z2 — »• Zitt), which is mainly a reflection of the 
increasing importance of virtual Higgs bosons in the Z2 three-body decays. We mention 
that for values of tan/? somewhat below the range where the decay Z2 Zih becomes 
kinematically allowed, contributions from all neutral Higgs bosons are important. 

The above considerations motivated us to begin a systematic exploration of how signals 
for supersymmetry may be altered if tan /3 indeed turns out to be very large. To facilitate this 
analysis, we have incorporated the above calculations into the computer program ISAJET 
7.32, so that realistic simulations of sparticle production and decay can be made for large 
tan p. 

Another important effect at large tan/5 is that tau Yukawa interactions can alter the 
mean polarization of the r's produced in chargino and neutralino decays. This, in turn, 
alters the energy distribution of the visible decay products of the r. The r polarization 
information is saved in ISAJET and used to dictate the energy distribution of the r decay 
products. 

The rest of this paper is organized as follows. In Sec. II, we describe aspects of our event 
generation and analysis program for Tevatron experiments, including a catalogue of some 
of the possible signals for supersymmetry at large tan/?. In Sec. Ill, we present numerical 
results of our generation of supersymmetric signals and SM backgrounds, and show the 
reach of the Tevatron MI and TeV33 in the parameter space of the mSUGRA model. In 
Sec. IV, we present a summary and conclusions from our work. Some lengthy three-body 
decay formulae are included in the Appendix. 

II. EVENT SIMULATION, SIGNATURES AND CUTS 

In several previous works , a variety of signal channels for the discovery of supersym- 
metry at the Tevatron were investigated, and plots were shown for the reach of the Tevatron 
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MI and TeV33 in the parameter space of the mSUGRA modeL The simulation of SUSY 
signal events was restricted to parameter space values of tan (3 = 2 and 10. The promising 
discovery channels that were investigated included the following: 

• multi-jet +I^rp events (veto hard, isolated leptons) (JOL), 

• events with a single isolated lepton plus jets (JIL), 

• events with two opposite sign isolated leptons plus jets (JOS), 

• events with two same sign isolated leptons plus jets +-^7^ (JSS), 

• events with three isolated leptons plus jets (J3L), 

• events with two isolated leptons +I$rp (no jets, clean) (COS), 

• events with three isolated leptons (no jets, clean) (C3L). 

In these samples, the number of leptons is exactly that indicated, so that these samples 
are non-overlapping. For Tevatron data samples on the order of 0.1 fb~^, the JOL signal 
generally gave the best reach for supersymmetry. It is the classic signature for detecting 
gluinos and squarks at hadron colliders. For larger data samples typical of those expected 
at the MI or TeV33, the C3L signal usually gave the best reach. In the present paper, we 
will extend these results to the large tan (3 region of mSUGRA parameter space; we will also 
look for new signatures which may be indicative of supersymmetry at large tan (3. 

By examining the branching fractions in Fig. |^, we expect in general at large tan (3 that 
there would be a reduction in supersymmetric events containing isolated e's or /I's. We 
also expect for large tan (3 and small mo a more conspicuous presence of isolated r leptons 
(defined by hadronic one- or three- charged prong jets as discussed below). For large tan/? 
and large mo, we expect an increased presence of tagged 6-jets (defined by displaced decay 
vertices or by identification of a muon inside of a jet). For these reasons, we have expanded 
the set of event topologies via which to search for SUSY to include, in addition: 

• multi-jet +Ifrp events which include at least one tagged 6-jet (JOLB), 

• multi-jet +Ifrp events which include at least one tagged r-jet (JOLT), 

• multi-jet +Ifrp events which include at least either a tagged 6-jet or a tagged r-jet 
(JOLBT), 

• opposite-sign isolated dilepton plus jet events where at least one of the isolated 
leptons is actually a tagged r-jet (JOST), 

• same-sign isolated dilepton plus jet events where at least one of the isolated 
leptons is actually a tagged r-jet (JSST), 

• isolated trilepton plus jet +Ifrp events where at least one of the isolated leptons is 
actually a tagged r-jet (J3LT), 

• clean opposite-sign isolated dilepton events where at least one of the isolated 
leptons is actually a tagged r-jet (COST), 
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• clean isolated trilepton H--^^. events where at least one of the isolated leptons is actually 
a tagged r-jet (C3LT). 

We note that some of these event samples are no longer non-overlapping; for instance, the 
JOLB sample is a subset of the canonical -^j. (JOL) sample. In the tau samples, the lepton 
multiplicity is again exactly that indicated, except that at least one of the leptons is required 
to be identified as a r. 

To model the experimental conditions at the Tevatron, we use the toy calorimeter sim- 
ulation package ISAPLT. We simulate calorimetry covering —4 < rj < A with cell size 
Arj X Acf) = 0.1 X 0.0875. We take the hadronic (electromagnetic) energy resolution to be 
70%/a/E' {15%/^/E). Jets are defined as hadronic clusters with Et > 15 GeV within a 
cone with AR = = 0.7. We require that |?7j| < 3.5. Muons and electrons are 

classified as isolated if they have px > 5 GeV, \ri{i)\ < 2.5, and the visible activity within 
a cone of R = 0.3 about the lepton direction is less than max{ET{i)/4:, 2 GeV). For tagged 
6-jets, we require a jet (using the above jet requirement) to have in addition |?7j| < 2 and to 
contain a 6-hadron. Then the jet is identified as a 6-jet with a 50% efficiency. To identify 
a r-jet, we require a jet with just 1 or 3 charged prongs with px > I GeV within 10° of 
the jet axis, and no other charged prongs within 30° of the jet axis. The invariant mass of 
the 3 prong jets must be less than m^, and the net charge of the 3 prongs should be ±1. 
QCD jets with px = 15 (> 50) GeV are mis-identified as r jets with a probability |jT6| of 
0.5% (0.1%), with a linear interpolation in between. In our analysis, we neglect multiple 
scattering effects, non-physics backgrounds from photon or jet misidentification, and make 
no attempt to explicitly simulate any particular detector. 

We incorporate in our analysis the following trigger conditions: 

1. one isolated lepton with pt(^) > 15 GeV and I^j, > 15 GeV, 

2. > 35 GeV, 

3. two isolated leptons each with Et > 10 GeV and I^j. > 10 GeV, 

4. one isolated lepton with Et > 10 GeV plus at least one jet plus I^j. > 15 GeV, 

5. at least four jets per event, each with Et > 15 GeV. 

Thus, every signal or background event must satisfy at least one of the above conditions. 

We have generated the following physics background processes using ISAJET: ti pro- 
duction, ly+jets, Z+jets, WW, WZ and ZZ production and QCD (mainly from bb and cc 
production). Each background subprocess was generated with subprocess final state parti- 
cles in pT bins of 25 - 50 GeV, 50 - 100 GeV, 100 - 200 GeV, 200 - 400 GeV and 400 - 600 
GeV. 

III. THE REACH OF THE FERMILAB TEVATRON FOR MSUGRA 

We present our main results for the reach of the Tevatron for mSUGRA at large tan /3 
in the tjiq vs. mi/2 parameter space plane for Aq = and for tan [3 = 2, 20, 35 and 45. Our 
results are shown for > only. For small tan/5 ~ 2, the /i < results differ substantially 
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from the fi > results, and are shown in Ref. 0]. As tan/5 increases, the positive and 
negative fi results become increasingly indistinguishable. 

In Fig. ^ we show for orientation contours of constant rrig and rrig in the mg vs. mi/2 
plane. The bricked regions are excluded by either lack of appropriate electroweak symmetry 
breaking, or due to the fi or Wi being the LSP instead of the Zi. The gray regions are 
excluded by previous experimental sparticle searches, and the excluded region 0] is domi- 
nantly formed by the LEP2 bound that > 80 GeV [^]. The most noticeable feature 
of Fig. ^ is that the theoretically excluded region increases significantly as tan f3 increases. 
In the low rriQ region, this is due to the decrease in fi mass, making it become the LSP. 
The contours of rrig and rrig on the other hand are relatively constant and change little with 
tan/3. The region to the left of the dotted lines denotes where the decay modes Wi ^ fiv 
and Z2 —>■ TiT become accessible. 

As in our previous analysis of signals at low tan/3 [^], for channels involving jets, we 
require of all signals, 

• jet multiplicity, rijets - n^-jets > 2, 

• > 40 GeV, and 

• EtUi), Et{j2) > and > E^, 

where the parameter is taken to be E!^ = 15,40,60,80,100,120,140,160 GeV. This 
requirement serves to give some optimization of cuts for different masses of SUSY particles. 

We generate signal events for each point on a 25 GeV x 25 GeV grid in the mo — mi/2 
plane. For an observable signal, we require at least 5 signal events after all cuts (including 
those detailed below) are imposed, with N signal exceeding 5^Nbackground- Any signal is 
considered observable if it meets the observability criteria for at least one of the values of 
E^. In addition, we require the ratio of signal/background to exceed 0.2 for all luminosities. 

A. Reach via the JOL channel 

As in Ref. 0], for multijet+^j, events (JOL), we require in addition to the above, 

• transverse sphericity St > 0.2, and 

• A(l){$^, Etj) > 30". 

In Fig. ^, we show the Tevatron reach via the JOL channel. Black squares denote param- 
eter space points accessible to Tevatron experiments with 0.1 fb~^ of integrated luminosity 
(approximately the Run I data sample); points denoted by gray squares are accessible with 2 
fb~^ while those with open squares are accessible with 25 fb^^. Points denoted by x are not 
visible at any of the luminosity upgrade options considered. In frame a) , no black squares are 
visible; regions normally accessible to Tevatron experiments with just 0.1 fb~^ of integrated 
luminosity have been excluded by the negative results of LEP2 searches for charginos. This 
is strictly valid only within the model framework, and should not be regarded as a direct 
bound on rrig. Regardless of the LEP2 bounds, Tevatron experiments should directly probe 
this region via the independent search for strongly interacting sparticles. Note that even 
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within the mSUGRA framework, for /i < and tan (3 = 2, where fTi^^ is considerable heavier 
for the same mi/2 values, there still exist parameter space points accessible with only 0.1 
fb~^ A significant number of gray squares appear in frame a), denoting regions with 
rrig ~ 400 GeV that can be probed at the MI. As tan f] increases, the theoretically excluded 
region absorbs some of these points at low mo, while some of the high mo points become 
inaccessible. In the latter case, much of the signal actually comes from WiWi and W1Z2 
production, and these particles decay decreasingly into jetty final states, so the JOL signal 
diminishes. Finally, for very large tan j3 = 45, none of the parameter space in this channel 
is open to MI searches. For TeV33, we see that mi/2 ~ 175 GeV (m^ ~ 475 GeV) can be 
probed in all of the frames a)-d) as long as mo is not much larger. The largest reach occurs 
when E!j. attains its largest value of i?^ = 160 GeV. 

B. Reach via the JOLB channel 

In Fig. 1^, we show the reach in the ^j^+jets channel, where in addition we require at least 
one tagged 6-jet (JOLB). Comparing with Fig. ^, we see that the requirement of a tagged 
6-jet considerably reduces the reach of the MI. Furthermore, the parameter space points 
with mi/2 = 175 GeV are no longer accessible to TeV33. In other words, a higher value 
is more efficient in maximizing signal-to-background for large mi/2 than requiring an extra 
6-jet. However, for large mo and mi/2 ~ 125 — 150 GeV, the extra 6-tag does somewhat 
increase the reach of TeV33 for SUSY. Comparison of Fig. ^ and § shows three additional 
points accessible in frame a), two in frame b), and one in frame d). We have also tried to 
extend the parameter space reach by requiring an identified r-jet (JOLT) or either a r or 6 
jet (JOLBT) along with jets. In both of these cases, no additional reach was achieved 
beyond the results of Figs | and ^. 

C. Reach via the JOS and JSS channels 

The reach of Tevatron upgrades on the JOS channel is presented in Fig. ^ We require, 
in addition to the conditions at the beginning of this Section, 

• events with exactly two opposite sign isolated leptons (e and /i), with -Et(^i) > 10 GeV 
and a veto of r-jets. 

At the Tevatron at low tan/?, signals in this channel mainly come from W1Z2 production, 
where Z2 decays leptonically, and Wi decays hadronically, while top production is a major 
source of SM background. There is significant reach by the Tevatron MI and TeV33 in this 
channel at low tan/5, as seen in frame a). As tan/? increases, the Z2 leptonic branching 
fraction decreases (see Fig. so that the MI has no reach in this channel for tan/5 > 20. 
The reach of TeV33 is severely limited in this channel at high tan f3 as well. 

We have also examined the reach of the MI and TeV33 for same-sign dileptons (JSS 
channel), where we require in addition 

• events to contain exactly two same sign isolated leptons, again with Exih) > 10 GeV 
and a veto of r-jets. 
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The reach of Tevatron upgrades in this channel for mSUGRA is not very promising. The 
signal should result mainly from gg and gq production mechanisms, but these have only 
small cross sections for parameter space points beyond the reach of LEP2. We found almost 
no reach for mSUGRA in this channel beyond the LEP2 bounds for any values of tan (3. 

We have also studied the Tevatron reach in the dilepton plus jets channels where we 
required in addition that at least one of the leptons be a tagged r-jet: the JOST and JSST 
channels. In each of these small increase in reach was obtained for large values of 

tan j3 and low mo beyond the corresponding "tau-less" channels. Most of this additional 
region can also be probed via the J3L channel discussed below, so we do not show these 
results here. 

D. Reach via the J3L channel 

For small values of tan /?, the J3L channel considerably increases the region of mSUGRA 
parameters beyond what can be probed via the channel at a high luminosity Tevatron. 
In addition to the generic cuts for all the signals involving jets, we require the following 
analysis cuts for the J3L channel: 

• events containing exactly three isolated leptons with ET{ii) > 10 GeV and a veto of 
r-jets, plus 

• we veto events with \M - Mz\ < 8 GeV. 

The reach in the J3L channel after all cuts are imposed is shown in Fig. |^. Since the signal 
almost always involves a leptonically decaying Z2, it is not surprising to see that the large 
reach at low tan/? is gradually diminished until there is almost no reach for tan/3 ~ 45. 

We have also examined the Tevatron reach in the trilepton plus jets channels where we 
required in addition that at least one of the leptons be a tagged r-jet: the J3LT channel. As 
before, only a slight additional reach was obtained at large tan /? and low mo beyond what 
could be probed via the "tau-less" J3L channel. Here, and in the jetty dilepton channels 
mentioned above, this is presumably because secondary leptons from tau decay tend to be 
soft, and fail to satisfy the acceptance requirements. Again, we do not show these results 
here. 

E. Reach via the C3L and C3LT channels 

For small tan/? ~ 2, and a large enough integrated luminosity, the maximum reach of 
the Tevatron was often achieved via the clean trilepton channel from W1Z2 ^ 3£ + Ifj,. For 
the C3L signal, following our earlier analysis we implement the following cuts: 

• we require 3 zso/atec? leptons (e and /i) within \ri(\ < 2.5 in each event, with Exiii) > 20 
GeV, ET{i2) > 15 GeV, and Erih) > 10 GeV, 

• we require I^j, > 25 GeV, 

• we require that the invariant mass of any opposite-sign, same flavor dilepton pair not 
reconstruct the Z mass, i.e. we require that |m(££) — Mz\ > 8 GeV, 
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• we finally require the events to be clean, i.e. we veto events with jets. 

Our calculated background in this channel is 0.2 fb. 

In Fig. 1^, we show the reach in the C3L channel for the four cases of tan f3. In frame 
a), we see at low tan/5 that indeed there is no reach beyond the current LEP2 bound in the 
C3L channel for 0.1 fb~^. For the MI integrated luminosity, however, there is considerable 
reach to values of mi/2 ~ 225 GeV, and for TeV33, the reach extends to mi/2 ~ 250 GeV, 
corresponding to rrig ~ 700 GeV! As tan/? increases, the branching fraction for a leptonic 
decay of Z2 and Wi decrease. In frame b), in fact, we find no reach for SUSY via the C3L 
channel for MI and considerably reduced reach for TeV33, except at large mo. For smaller 
values of mg a complicated interference between various amplitudes reduces the leptonic 
decay width of Z2. As tan/? increases even further to 35 and 45 as in frames c) and d), the 
C3L reach is wiped out at low mo- Some reach remains at large mo in frame c), where the 
branching fraction BF{Z2 iiZi) ~ BF{Z It). In frame c?), most of this region also 
becomes inaccessible because of the increased importance of (virtual) Higgs boson mediated 
decays of Z2 which lead to a strong enhancement of its decay to hhZi . 

We have also examined the reach for clean trileptons, where one of the leptons is actually 
an identified r-jet (C3LT). In this case, we relax the additional px requirements on the 
leptons. This increases the chance of detecting the softer secondary leptons from the decay 
of tau(s). Trigger 4 presumably plays an important role for this class of events. The reach 
via this channel is shown in Fig. ^ In frames fe), c) and (i), significant additional reach is 
gained in the low mo regions, beyond that shown in any of the previous figures! Notice that 
the region where the signal is observable is where chargino and neutralino decays to real fi 
are accessible (see Fig. The reach in the C3LT channel effectively extends the reach of 
TeV33 to mi/2 ~ 250 GeV for at least some value of mo for all the values of large tan/? 
considered. We remark that the gain in reach via channels involving taus is limited because 
we require the presence of additional hard leptons (e or /i), jets or I^j, in order to be able 
to trigger on the event. Because secondary leptons from the decay of a tau tend to be soft, 
the development of an efficient r trigger may significantly enhance the reach when tan /? is 
large. 

F. Reach via the COS and COST channels 

In our previous studies [Q we had already noted that for small values of tan /?, a study 
of the clean opposite sign dilepton channel (COS) would allow a confirmation of the signal 
in the C3L channel for a large range of mSUGRA parameters. For the COS channel, we 
require 

• exactly two isolated OS (either e or /i ) leptons in each event, with Et{C.i) > 10 GeV 
and -Et(^2) > 7 GeV, and \r]{£)\ < 2.5. In addition, we require no jets, which effectively 
reduces most of the ti background. 

• We require I^j, > 25 GeV to remove backgrounds from Drell-Yan dilepton production, 
and also the bulk of the background from 7*, Z ^ rf decay. 

• We require (p{£i) < 150°, to further reduce 7*, Z ^ rf background. 
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We require the Z mass cut: invariant mass of any opposite-sign, same flavor dilepton 



pair not reconstruct the Z mass, i.e. m{££) — Mz 
B = \0^\ + \pT{h)\ + \pt{^^)\ < 100 GeV. 



> 8 GeV. Finally, we require 



Our calculated background in this case is 64 fb. 

We have checked that while there is an observable signal at the MI (TeV33) for mi/2 ~ 
150 (175) GeV, and if uiq < 100 GeV, there is no observable signal for any of the allowed 
regions of the plane if tan/3 > 20. We have also examined this channel by requiring in 
addition that at least one of the leptons be an identified r-jet (COST). In this case, no reach 
for mSUGRA was found for any of the tan /3 values considered. We therefore do not show 
these figures. 



IV. SUMMARY AND CONCLUSIONS 

To summarize the reach of Tevatron upgrades for large and small tan /5, we show in 
Fig. ^ the SUSY reach via all of the channels that were examined, for both the upgrade 
options of the Tevatron. Thus, if a parameter space point is accessible via any channel, 
we place an appropriate box, corresponding to the integrated luminosity that is required. 
The cumulative reach shown in the figure is completely established with just four channels: 
JOL, JOLB, C3L and C3LT. For some points, the signal may be observable in more than one 
of these or other channels studied in this paper. It is possible that some additional reach 
may be gained by combining several channels to gain a net "5cr" signal, even though the 
significance in each of these channels is somewhat smaller. We do not consider this added 
detail here. 



We see from Fig. |10| that as tan (3 increases, the SUSY reach of Tevatron upgrades is 
significantly reduced. For the MI option, there is no reach beyond current LEP2 bounds 
that can be established at tan/3 = 45. The TeV33 option has some reach in all frames, but 
clearly a much reduced reach for large tan j3. In particular, there are parameter regions just 
beyond the current LEP2 bounds for which there will be no observable signal even with the 
luminosity of TeV33. The reduction of the reach is mostly due to the depletion of leptonic 
signals, especially the clean three lepton signal, in the region of large tan (3. Note that the 
branching ratio for Wi and Z2 to decay into electrons and muons plus missing particles is 
actually quite large if charginos and neutralinos dominantly decay into real or virtual fi. 
However, the secondary leptons produced in subsequent r decays are usually too soft to 
pass our trigger criteria or acceptance cuts. It might be worthwhile to investigate whether 
these cuts can be lowered without introducing unacceptably large backgrounds {e.g. from 
heavy flavors, where the lepton happens to be isolated and the jet is lost, or from jets faking 
leptons) or via a development of a special trilepton trigger. 

Modes with identified (hadronically decaying) taus could only partly compensate this 
loss of reach in the leptonic channels. Again the problem seems to be that the hadronic 
decay products of the r leptons are frequently too soft to pass the cut Et{t — jet) > 15 
GeV. It might be worthwhile to study if this cut can be lowered, e.g. by focussing only on 
one-prong r decays, for which QCD backgrounds are much smaller than in the three-prong 
channel. In addition, the triggers adopted in our study are not very efficient for events 
with rather soft leptons plus r— jets, as in our C3LT sample. We therefore believe that the 
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reach of future Tevatron runs could be extended significantly in the region of large tan (3 if 
it is possible to devise strategies to reliably identify, and perhaps even trigger on taus with 
visible pt smaller than 15 GeV. We remark, however, that even without such developments, 
experiments at the LHC will probe the entire parameter plane shown at least via the 
channel. 
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In this Appendix we give analytical expressions for those three-body partial widths that 
are sensitive to 6 or r Yukawa couplings and/or to — fn mixing (/ = 6, r). We first list 
the relevant couplings, and then give results for Zi — > Zjff, Wi — > rz/,-Zj, Zi Wjff, and 
~9 ^ btW,. 

Many of the couplings and kinematic functions that enter our computations have been 
defined in our earlier papers. Instead of rewriting these lengthy definitions again, we provide 
the reader with references to the papers from which these couplings are used. In these studies, 
the two charginos were denoted by W_ and W+ instead of Wi and W2, respectively. Also, 
the lighter (heavier) neutral CP even Higgs scalar was denoted by Hi (Hh) rather than by h 
(H), while the CP odd pseudoscalar was denoted by Hp rather than A. The corresponding 
couplings are characterized by superscripts I, h and p. To facilitate the use of these couplings 
from the earlier literature, we use this older notation to denote the charginos and neutral 
Higgs bosons in the formulae listed in this Appendix. 



The couplings of electroweak neutralinos and charginos to a fermion and a sfermion are 
affected by mixing between SU{2) doublet (L— type) and singlet (-R— type) sfermions. We 
write the sfermion mass eigenstates as: 



where /i denotes the lighter eigenstate. Since there is no L — i? mixing in the sneutrino 
sector, some couplings remain unaffected. We list these for completeness, using the notation 
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Appendix: Sparticle Decay Widths for Large tan/3 



1. Couplings 



= cosOffi - sine ffn; 
/a = sinOffL + cosO ffn, 



(1) 



of Refs. 0] and [|T|]: 




(2a) 
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^^_ = -^sm7R; (2b) 
= sin (2c) 
= -/.cos7L. (2d) 

Here, g and g' are the SU{2) and f/(l)y gauge couplings, and // the Yukawa couphngs of 
fermion /. The corresponding couphngs of the heavier chargino mass eigenstate Wj^ can be 
obtained by the substitutions 



W_^W+: cos'jL,R^ -Ox,ySm'yL,R; sm'yL,R^ Ox^yCOS'yL,R- (3) 

In the calculation of the partial widths, we will ignore terms oc mb,mr when doing the 
Dirac traces. It then becomes convenient to write the matrix elements in terms of couplings 
to fermions with fixed chirality. In the following we denote all left-handed couplings with 
the symbol a, and right-handed couplings with f3. The chargino couplings to the lighter 
third generation squark mass eigenstates can be written as: 

= -g sin cos 9t + ft cos 7/? sin Ot] (4a) 

~ _ = - /b cos 7l cos ; (4b) 

a ~_ = -g sin 7i cos 6b + fb cos 7^ sin Ob, (4c) 

4_ =-/tCOS7«cos^,. (4d) 

The corresponding couplings to third generation sleptons can be obtained by the substitu- 
tions: 

q^l: e^^O; Ob ^ 6,; f^^O; h^fr- (5) 

Similarly, the couplings to the heavier sfermion mass eigenstates /2 can be obtained by 
substituting: 

fi—>-f2'- cosOf sinOf] sm6f^—cos6f. (6) 

Finally, the couplings of the heavier chargino state can again be computed using Eq.(|^). 
The couplings of neutralinos to b and r (s)fermions can be written as: 

al^ = ALcos9f-ffV^hm9f- (7a) 

4] = ff^2^ cos ^/ + bL sindf, (7b) 

where A~ , B~ are as in Ref. |]19l. The couplings to fermions with weak isospin 1^ = +1/2 
can be computed from eqs. (|^) by inserting the corresponding unmixed couplings; in addition, 
one has to replace the component of the neutralino eigenvector by The couplings 
to heavier sfermion eigenstates can again be obtained by applying Eq.(|D. 
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Finally, we introduce the charged Higgs - chargino - neutralino couplings 



4 ) 



where i is the neutralino index; the can be found in Ref. |20|. 



(0 



(3g = -sinPA^i'e, 



(8a) 
(8b) 



2. Z, ^ ZJf Decays 



We are now in a position to present our results for the partial widths for decays involving 
third generation fermions. We begin with the decay of a neutralino into a lighter neutralino 
and a bb or r+r^ pair. This decay can proceed through the exchange of the two sfermion 
mass eigenstates fi^2, through the exchange of a Z boson, or through the exchange of one of 
the three neutral Higgs bosons of the MSSM. The partial width can therefore be written as 



1 



r(Z,- ^ Ziff) = -NJf)-^ 



(9) 



where the color factor Nc{f) = 3 (1) for f = b (r). Recall that we set m/ = when 
evaluating Dirac traces. As a result, the Higgs and Z exchange diagrams do not interfere 
with each other.Q 

The pure sfermion exchange contribution is given by 



/2 /l,2' 



where 



LL 



fk 

LR 



(A; = 1,2); 



-'/lp/2 _|_ p/lp/2 



/if/2 



/if/2 



(10) 

(11a) 
(lib) 



Here, the subscripts L and R refer to the chirality of the SM fermion coupling to the heavier 
neutralino Zj. The quantities appearing in Eq.(|rT|) are: 



LL 



4 a 



a 



,2 / - X 2' 
fk\ _j_ j pfk 



(12a) 



^This would be a very bad approximation for Zj Zitt decays. However, if these decays are 
allowed, Zj has numerous 2-body decay modes into real gauge and Higgs bosons and lighter 
neutralinos and charginos. The branching ratios for neutralino 3-body decays into top quarks are 
therefore always negligibly small. Analogous remarks apply to Wj — > Zitb decays. 
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-pfk 

^ RR 



4/3 



a- 



LR 



«f/i«t/^t^("^^,' "^z,); 



Zo Zi Zi Zi 



Zi Zi Zi Zi 



3/1 af2 



-8«|/^|"f/t^(^z- '^/i' "^/.' '^z,); 



(12b) 
(12c) 

(12d) 

(12e) 
(12f) 
(12g) 



The kinematic functions ip, 0, and Y are given in Ref. [^0, and 6i is (1) if the sign of 
the i*^ eigenvalue of the neutrahno mass matrix is positive (negative) [0. The functions 
ip and 0, which depend on two sfermion masses are generahzations of the functions and 
(j) which depend on just one sfermion mass: to define ip, we simply split the squared factor 
where the stop mass occurs in Eq. (A6a) of Ref. |2l|, into two such factors, with each one 
containing a different sfermion mass. Similarly, (p is generalized from (p by writing mj_^ in 
the first factor outside the square parenthesis in Eq. (A6b) of Ref. |21], and mj^ inside the 

square parenthesis. In other words, when the two sfermions /i and /2 have the same mass, 
ip = ip and (p = (p. 

For completeness, we also give the squared Z exchange contribution, which is not affected 
by sfermion mixing: 



2 



dE- 



m- 



M| - 2Em7 



E 



Zi Zj ^ ' Zi Zj 



E^ + m% + ^(E'^ - m 



^Note that the third line in Eq.(A6h) of that paper should come with a positive overall sign. 
Furthermore, the last term in the first denominator in Eq.(A6a) should be m?, rather than m^. Of 
course, nit is repaced by the appropriate fermion mass in the definition of these functions. Finally, 
although the number of arguments appearing in the Y function are different from that in Ref. |21|, 
the correspondence is obvious. 
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2mi 



(13) 



Here, e is the QED coupling, Wij is the ZZiZj couphng given in Ref. ||2^, and «/, /?/ are 



the left- and right-handed Zff couplings in the notation of Ref. [^. Finally, the upper 
integration limit is given by 



m~ + m~ — 4m f 
2m^ 



(14) 



and 



\ 



m~ + m~ — 2Em'r 



(15) 



The ]9ure scalar Higgs exchange contribution can also be written as a single integral: 

2 „ , 



2n' 



gnif 



M]v cos/5 



2 2 

m~ + m~ 



Z^ 



2m;? — 2m 



sma 



m~ + m^ 



2m7; E — m?„ 
Zj Zj 



+ 



z. 



'm; 



COS a 



m~ + m~ — 2m;y — m%r 

Zi Zj Hh 



z^ 



(16) 



Here, Xlj^ are the couplings of the light and heavy neutral scalar Higgs boson to two neu- 
tralinos and a is the angle describing mixing in the scalar Higgs sector as defined in Ref. 
20| , and m^^ ^ are the masses of the two Higgs bosons. The upper integration limit is again 



given by Eq.(0) 

The squared pseudoscalar Higgs exchange contribution can be cast in a quite similar form: 



r 



Hr, 



27t' 



En 



gmfta.n(3 



w 



1 2 



m; 



dEBfjE'^ - m| 



m~ + m~ 



2m~ E — 2m 



z^ 



m~ + m~ — 2m ;j E — m'ir 



(17) 



The couplings Xf can again be found in Ref. pO 



We now turn to the various interference terms listed in Eq.(^. The Z—sfermion inter- 
ference contributions can be written as 



r 



Zf 



r - + r - 



(18) 



17 



with 

r 



jfk afk 



2m 7 



( V 



ds 



m~ — s — m^f 



Am-r 

m^^ [Eq + Q') - 
■ m~ {Eq - Q') - 

Here we have introduced the quantities 



2 I 2 2 

Jk Zi 



s + m~ — m~ 
2m7i 



Q = \/El - s, 



(19) 



(20) 



and 



Q' = Q\i 



Amj 



The real couphng Wmn is defined to be 



-I 



with Wmn given in Ref. [p^ ]. 

Finally, the Higgs-sfermion interference contributions can be written as 



(21) 



(22) 



where Hi, Hh and Hp again denote the light scalar, heavy scalar, and pseudoscalar Higgs 
boson, respectively. The separate contributions in Eq. (E^ are given by: 



27r^ gmf sin a 



Jh {m^ , m/^ , mm , , 6^^ + Oj ) 



Hhfk 



271^ gmf cos a 
m^ cos/5 



(-1)^- 



(m^^ , mj^ , m^,, , m^^ ,ei + 9j); 



lip Ik 



2tt^ gmftanP / j 



m^^ Mw 



Jh (m^^ , mj^ , m^^ , m^^ , 1 + 6*^ + 6* 



(23a) 



(23b) 



(23c) 
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The function Jh is defined as 

f j m m~ J (J^g 

Juimj ,m7,mH,mj ,9) = ^ (24) 

^■J ' Jim) s - mjj 

'l smj- - m){m^^ + m^^) + (-l)^m^^m^^ (s - 2m)) {Eq + Q') - /i^" 

where /z^, i^g, Q and Q' have been defined in Eq. (pUD . 



3. Wj — > ZiTVr Decays 

These decays proceed via the exchange of a boson, a charged or neutral third gener- 
ation slepton, or a charged Higgs boson. The partial widths can thus be written as 



1 1 



(25) 



The Higgs and W exchange contributions do not interfere, since we neglected terms oc rrj-T- 
when doing the Dirac algebra. 

The squared W exchange contribution is given by 



2 rE^ 



V = / dE- 



^2 _ 



-3 



X] 



+ 



1^^ 



3(m| +m| Jm^ 



E-2m~ f2E^ + m| 



m,7, ( m~ + m~ — 2Em 



(26) 



Here X] and F/ are the WWjZi couplings as defined in Ref. |]T8[, and the upper integration 



limit ii^max is given by Eq.(|T^) with 



myT, and m/ — >• 0. 



The squared sneutrino exchange contribution is given by 



r. = 2(I|)^ 



(27) 



The couplings appearing in Eq.(P7D have been defined in eqs.(0), and the kinematical func- 
tion ip is given in Ref. . 



The pure scalar tau exchange terms can be written as 

r- = r- + r- + r- 



T2 I ^ TIT2 1 



(28) 
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where 



2 a 



w., 



a 



Zi Zi Zi Zi 



(29a) 
(29b) 



The couphngs appearing in eqs.(p9D have been defined in eqs.(||)-(^, and the functions ip, ip 
are as defined above. 

The squared charged Higgs boson exchange contribution is 



L H = m 



' gmrtan(3\' 
, Mw I 



(30) 



E 



a 



+ A 



Wi 



m ~ + m~ — 2£'m 



Here, -Emax is the same as in Eg. (P^) , the H^W^Z-i couphngs have been defined in Eqs.(||), 
and 6j (= 6'_ or 6*+ in the notation of Ref. [|18[) = (1) if the corresponding eigenvalue of 
the chargino mass matrix is positive (negative). 

The W—sneutrino interference contribution is not affected by — mixing and con- 
tributions (X fr'i it can be written as 



Tw, = -AV2g'i-lf'-'^A^~A''~ 



(31) 



where we have introduced the functions 

^2 



/i(mg„mp m^) 



TT 



2m^ Jo 



ds 



w 



f 



mir — m~ — m 



f 



w 



4m 



/2(m^,mp m^) 



JO 



log + - 



w 
dn 

"mi:'"^' (Eq - g) - 



(32a) 
(32b) 



m^, mg: ^ and mj^ 



The quantities /x^, Eq and Q are as in Eq . (120|) , with 

The same functions also appear in the W— scalar tau interference contributions: 

where 



m 



f 



(33) 



(Xj + F;) Ji(mg.^,,m^,,m^j - (Xj - F/) /2(mg.^., m^,„ m^j] . (34) 
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The couplings Xj, YJ can be found in Ref. [|T^]; the remaining couphngs appearing in Eq.(p4| 
have been introduced in eqs.(|^)-(0). 

The sneutrino- scalar tau interference terms can be written as: 



Tuf, + r 



UT2 } 



where 

r 



(35) 



(36) 



The functions Y, cf) have aheady been defined. 

The charged Higgs-sneutrino interference term is given by: 



where we have introduced the function 

^2 



1 r 9(0 



(m~ —m~ ) 
w ■ z- 



ds 



s — mjj 1 2 



(37) 



(38) 



4m 



w, f 



1q ^w, (^q + Q)-f^^ 

^w. {Eq -Q)- fi"^ 



the quantities /i^, i^g and Q are as in Eq. (pOD , with - The charged Higgs 

couphngs appearing in the integrand of Eq. ( |55D have been defined in eqs. (|) . 

The same function also appears in the charged Higgs-scalar tau interference contribu- 
tions: 



Ht 



r^f, + r 



Hf2-, 



where 



F o ^fc ^f, ^m^tan/3 



(39) 



(40) 



The partial widths for the analogous neutralino to chargino decays are given by crossing: 
T{Zi W;r-Ur) = T{Wf ^ Z,T-iyr){m^^ ^ m^). (41) 



Note that Zi can also decay into Wfr^Vr final states, with equal probability. However, 
these neutralino decays are usually not very important, since they are either phase space 
suppressed, or have to compete with 2-body decays of the heavy neutralinos. 
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4. g ^ With Decays 

These decays proceed through the exchange of any of the four stop and sbottom mass 
eigenstates; in the hmit db, fb ^ considered in the existing hterature |21], only one of 
the two sbottom eigenstates contributes here, since Br does not couple to charginos in this 
limit. Fortunately the general case does not introduce terms with new Dirac structure in 
the matrix elements; the necessary phase space integrals can therefore be extracted from the 



Appendix of Ref. 



We begin by defining eight kinematical functions: 

dEtPtEt (m? + - 2Etmg - 



G2{mg,m-^,m^) 



Gzijn-g.mi.m^) 



Gfi{mg,mi^,mi^,m^) 



m 



w 



(m'^ + mf — 2Etmg — (m'^ + mf — 2Etmg 
dE-bEl\^'\ml + ml- 2E-bmg, m\, m?) 



(42a) 



+ — — 2EbVrig — m 



w 



(m^ + ml — 2Ebmg — m?^ (^m| + ml — 2Ebmg 
{-ifwAmgm^mt J dEiElX^I'^iml + ml - 2E-bm~g, m|, ml 



(42b) 



(m^ + ml — 2Ebmg — m?^ (^m? + ml — 2E] 

dEt 



-1) « ^mgm^ 



m| + mf — 2Etmg — 



£"5 (max) — i?^(min) 



m? + ml — 2Etmg — 



2m?, 



mj + ml- 2Etm~g - vr?- 

2 ^ ; 



* m| + — 2Etmg — m'j 



(42c) 

(42d) 
^logX ; 

(42e) 



dEt 



2 J m| + m^ — 2Etmg — m 



I mg (^m? + mf — 2Etm 



9 



ml — ml / „ 

^ [2Etmg - ml - 



m- 



m-g 



+2 (2Etmg -ml - ml 



logX 

[E-b{m^x)-E-b{mm)]}; (42f) 



•1) 



dEt 



m^ + mj — 2Etmg — m^ 



2 [E5(max) - £;;(min)] ^— ^ logX ^ ; 



(42g) 



'lY'^mtmg 
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(^i + - 2Etmg - m^) [^'^(max) - Eiijmxi)] ^^^^^ 
{rnjj + ml — 2Etmg — "m^^) + — 2Etmg — rnf^^ 

Here, Og = (1) if a positive (negative) gluino mass parameter is chosen, and 6^ (= 6_ or 6'+ 
in the notation of Ref. [1^) is the corresponding quantity for the chargino mass eigenstate. 
Further, we have introduced i?5(min, max) which are given by 



m? + mf — 2mgEt + ml — m'^~){mg — Et) =]= pt\^/'^{m'^ + m^ — 2mgEt, ml, m^ 



w 



w' 



2 {m}j + m} — 2Etmg 



Also, 



and 



(43) 

(44a) 
(44b) 

Finally, the limits of integration over Et in 

eqs.(^) are from mt to (^m? + ml — {m^ + m;,)^^ /2mg, while the integration over Ei in 

eqs.(^^,c) goes from m^, to m? — (^m^ + /2m^. 

The partial widths for the processes under consideration can be written as 



X 



m? + 2Ei{m.ax)mg — m| 
m? + 2i?5(min)m^ — m| 



r(^ thWi) 



{2nY 2mg 



-TX 



(45) 



k,l=l 



where Qs is the S'f/(3)c gauge coupling. Note that in the limit m^ —>■ the two sbottom 
exchange diagrams do not interfere with each other. The individual contributions in Eq. 
are given by: 



4j + fe/' 



Gi{mg,mi^,m^J 
-(-1)'= sin(2^JG8("^5, 



hh = -2 [a%a%+(3lf2] cos(2^,)G8(m^, m,-^, m^j; 



a 



G2{mg,mz,m^) - a'^i /3''-iG3{mg,mz,m^); (46c) 



cosOt sine^a^P^^ + sin6't cosOi^p^ai^^ ) GQ^m-g, mi^,mi^,m^^ 



cosOt cos6'6a^^a*i^ + sin6't sin6'b/5^^/5*i J G^im-g, m^-^, m^^, m^^,) 
+ ( cos6't cos6'fe/3^^a*i^ + sin6't smOba^^^p^^^ \ Gz{mg, mi^ , m^^ , m^) 



- ( cos6't sui9b(3^ (3~^ + sin6'i cos6'ba~ a~ ] G'j{m-g, mj^, m^^, 



(46a) 
(46b) 



(46d) 
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The couplings appearing in eqs.(|4^) are listed in eqs.(^)-(|^). The other stop-sbottom inter- 
ference contributions can be obtained from Eq. ([46 d| ) by substituting the appropriate coupling 
constants and squark masses; in addition, one has to apply the substitution rules to the 
factors in Eq.( [46d|) that depend on third generation squark mixing angles. Finally, we note 
that gluinos have the same partial widths for decays into tbW~ and ibW^ final states. 

These formulae have been incorporated into the event generator IS A JET 7.32 |Q. Finally, 
we remark that we have also updated the formula for r(^ ^ tiZi) that appears in Ref. 
to include ti — mixing effects. This has also been incorporated into ISAJET. 
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FIG. 1. Selected sparticle and Higgs boson masses versus tan/3 for the mSUGRA model for 
parameters a) (mo,mi/2)^o) = (150,150,0) GeV and b) {mo,mi/2, ^o) = (150,500,0) GeV, for 
both signs of the parameter n. We take mt = 170 GeV. 

FIG. 2. Chargino (Wi) and neutralino (Z2) branching fractions versus tan/?. In a) and 
b), we take the parameters (mo,mi/2i^o) = (150,150,0) GeV while in c) and d) we take 
(mo,mi/2,^o) = (150,500,0) GeV. In all frames, fi > and rrit = 170 GeV. 

FIG. 3. Gluino and squark mass contours in the mo vs. mi/2 parameter plane, for a) tan (3 = 2, 
b) tan/3 = 20, c) tan/3 = 35 and d) tan/3 = 45. In all frames, we take ^0 = 0, ;U > and 
mt = 170 GeV. The bricked regions are excluded by theoretical constraints, while the gray regions 
are excluded by LEP2 bounds on . 

FIG. 4. A plot of points accessible to Tevatron MI and TeV33 searches for mSUGRA via 
multijet events. A 5a signal above background is found for some value of Ej. for the MI for gray 
squares, while white squares are accessible only at TeV33. Points with a x symbol are inaccessible 
to MI and TeV33. 

FIG. 5. Same as Fig. ^, except we require in addition that at least one of the jets be an identified 
6-jet. 

FIG. 6. A plot of the reach of the Tevatron MI and TeV33 for mSUGRA via the JOS signal. 
FIG. 7. A plot of the reach of the Tevatron MI and TeV33 for mSUGRA via the J3L signal. 
FIG. 8. A plot of the reach of the Tevatron MI and TeV33 for mSUGRA via the C3L signal. 

FIG. 9. A plot of the reach of the Tevatron MI and TeV33 for mSUGRA via the C3LT signal. 

FIG. 10. A plot of the combined reach of the Tevatron MI and TeV33 for mSUGRA via all of 
the signal channels considered in this paper. 
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